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We study the equilibrium properties of an Ising frustrated 
lattice gas with a mean field replica approach. This model 
bridges usual Spin Glasses and a version of Frustrated Per- 
colation model, and has proven relevant to describe the glass 
transition. It shows a rich phase diagram which in a definite 
limit reduces to the known Sherrington-Kirkpatrick spin glass 
model. 



I. INTRODUCTION 

The Ising Spin Glass transition has been discovered 
to describe many seemingly different phenomena, and to 
model real systems much beyond what originally thought 
1^. The idea of introducing frustrated Hamiltonians to 
capture the essential physics of glasses, random media 
properties, evolutionary models, protein and RNA fold- 
ing, granular packing, dynamics of complex flow, and 
many others, has grown fertile in the last years |p|-pl|. 

Recently, in this panorama, an Ising Spin Glass like 
model, a general version of the frustrated lattice gas, has 
been introduced for its new interesting Monte Carlo dy- 
namical and equilibrium features relevant to the descrip- 
tion of the glass transition |^,|| , and with some variations 
has been applied in the context of phase transitions in 
granular packing 

In this paper we study the mean field equilibrium prop- 
erties of such a model, adopting standard replica formal- 
ism. The system we consider is characterized by the fol- 
lowing Hamiltonian: 
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Sij Si Sj 



where the lattice gas site variables rii — 0, 1 have an 
Ising internal degree of freedom Si — ±1 and {ij) denotes 
summation over all nearest-neighbor pairs of sites. The 
Eij = ±1 are quenched random variables, his a. magnetic 
field applied on the system and /i is a chemical potential 
for the site variables. Essentially, the model considers 
a lattice gas in a frustrated medium where the particles 
have an internal degree of freedom (given by its spin) that 
accounts, for example, for possible orientations of com- 
plex molecules in glass forming liquids. As stressed by 



Coniglio , these steric effects are greatly responsible for 
the geometric frustration appearing in glass forming sys- 
tems at low temperatures or high densities. As is detailed 
in ref. Q , this model offers a clear and intuitive picture of 
the mechanism leading to a glass transition, qualitatively 
reproducing the complex dynamical behaviour present in 
this regime. 

The presented Hamiltonian is a natural bridge be- 
tween Frustrated Percolation and standard Ising Spin 
Glasses (SG). Indeed, this two apparently different mod- 
els are obtained in two definite limits of its parameters. 
If /Lt ^ oo and J / fj, 0, for energetical reasons each 
site must be filled and the known Ising Spin Glass is 
obtained; on the contrary if J — s- cx) and ^/ J ^ 0, gener- 
ally even at T = the configuration with each site filled 
is impossible. In this last limit only site configurations 
which do not close "frustrated loops" , i.e. loops of filled 
sites whose spins are not satisfying all mutual interac- 
tions, are allowed. This corresponds to a site Frustrated 
Percolation in which clusters have a further weight fac- 
tor of 2^" [Nc is the number of clusters in the system). 
Moreover, introducing clusters a la Kasteleyn and For- 
tuin, this Hamiltonian may be described in terms of a 
site-bond correlated and frustrated percolation. 

With a simple transformation, = SiUi, this model 
may be changed into an Ising spin-1 Blume-Emcry- 
Griffiths model (BEG) jl^ in which only the bilinear 
coupling is affected by the quenched disorder £.y, while 
the spin biquadratic term has a coupling with opposite 
sign in relation to the original BEG model ||l^ (in no- 
tation of |13|, K/J — —1) and no disorder. In the case 
without frustration (e^j = 1), that has also been stud- 
ied in p3[ |, the order parameters are the diluted mag- 
netization TO — (Sn) and the particles density d = (n). 
In the /i — !■ cx) limit one recovers the Ising model, with 
d — 1. In mean field, at T = 0, d = to = 0(/i) if /i 7^ 
and d — 771 = 1/3 for /i = 0. This point with density 
< c? < 1 will become an interval when frustration is 
introduced. Moreover, expanding for small to we obtain 
the equation satisfied by the critical temperature 



— = 1 + exp f 1 - — 



(2) 



and one can see that this transition line is reentrant [|13| 
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effect that will also disappear when introducing frustra- 
tion. 

In the following sections we study the replica mean 
field theory of Hamiltonian with a gaussian dis- 
tributed coupling. The phase diagram of the model 
presents several interesting regions depending on the val- 
ues of T and /i/J. For highly negative values of /i/J, 
there is only a paramagnetic phase. Lowering the tem- 
perature at small negative values of ^/J, the system has 
a discontinuous transition to a spin glass phase, in which 
even at zero temperature the density is lower than one. 
Increasing /i/J, the spin glass transition becomes contin- 
uous, while the zero temperature density is still below 
one. This occurs up to a certain point where the den- 
sity becomes unity (at T — Q) signalling that we are ap- 
proaching the Sherrington-Kirkpatrick limit. The Parisi 
replica symmetry breaking solution seems to hold when- 
ever a spin glass transition is encountered. 



II. MEAN FIELD RESULTS 

We present in this section results obtained in the mean 
field approximation. The starting point is the calculation 
of the free energy / according with the replica trick [0 : 



(3f = - hm 



ln[Z"],,„. 



n^o nN 



(3) 



where [. . ]av. stands for the average over the disorder, 
which we suppose gaussian with zero mean and variance 
jViV. We obtain: 
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a<b 



- In Tr Q-P^-sf 



where the single site Hamiltonian is: 



(4) 



Hef f = -/3 J' labSVS'n' - J f ^ - 1 j ^ dan'' 

a<b ^ ^ a 

(5) 



The self consistent equations for the order parameters 
are given by the saddle points of / and read: 



da = (n'^) 



(6) 
(7) 



where the average is done using the effective Hamilto- 
nian. The overlap qab has a certain degree of dilution in 
respect to the parameter introduced in the SK model | |T^ ] 
and reduces to it in the limit d„ = I. 



A. Replica Symmetry 

To get a general qualitative picture of the phase dia- 
gram of the system, we first made a simple replica sym- 
metric (RS) assumption, that is, Qab = 9(1 — Sab) and 
da — d. The free energy then reads: 



Pfs = - \p^J^{q^ - d^) - ^pjd^ - In 2 

- J Vz \n[l + e" coshiPJ^z + f]h)] 



(8) 



where the gaussian measure is Vz 
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id-q)+f3{n-Jd) 



(9) 



The saddle point equations obeyed by the order param- 
eters are: 



d^ V 



and 



q= Vz 



sh{l3J^z + I3h) 
- cosh{f3Jy/qz + f3h) 



smh^{l3J^z + Ph) 
^ + cosh(/3J^z + l3h)Y 



(10) 



(11) 



The effects introduced by the magnetic field will in gen- 
eral not be considered here and in what follows we take, 
unless mentioned, h — 0. 

To characterize the system, in fig.|l| we present some 
representative curves for both q and d for several values 
of /i/J. For large values, the system approaches the SK 
limit (Tc — > J, d ^ 1 and q qsK)- For /i/J > —0.56, 
the system has a continuous transition {q ^ 0) at Tc 
satisfying 



— = 1 + exp 1 — 



(12) 



value to be compared with eq.(Q). Decreasing /i/J fur- 
ther, the transition line becomes first order (for —0.77 < 
/i/J < —0.56) as signaled by a jump in the order param- 
eter. At the transition line a partial freezing takes place 
(q < c? < 1), behaviour that has to be compared with 
other disordered models with discontinuous transitions 
like, for instance, the p-states Potts glass with p > 4 [T^ ] 
and the p-spin interaction model with p > 2 [|l6[ . We can 
also see from fig.0 that at low temperature q approaches 
d, while the actual value they assume at T = 0, where 
q — d (the system being fully frozen), depends on the 
chemical potential as can be seen in fig.0. The point 
where the transition changes behavior (/x/J ~ —0.56) 
turns out to be a tricritical one when including a non 
zero magnetic field. The system is a simple paramagnet 
{q = 0) below this region. This information is summa- 
rized in the phase diagram T versus /i presented in figj^. 
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The reentrant phase found in the case without frustration 
1^ is replaced here by these various regions. 

It is interesting to study the T = hmit of this model. 
For /i above the point 



J 



1 
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/27r 



0.6 



(13) 



it is possible to see that d — 1 and C = f3J{d — 
q) — •y/2/7r, results known to be characteristic of the 
Sherrington-Kirkpatrick SG. Something changes below 
H*, eq.(0), where the saddle point equations give: 
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C = 



■ exp 



hc + d- 
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(14) 
(15) 



These equations imply d < 1 at T = 0, as it should be 
in the true J ~f oo limit. On the other hand, below 
/ic — —0.77, just the paramagnetic solution q = d ^ 
is present, the transition being discontinuous (see fig.^, 
in accordance with the above phase diagram. The main 
novelties present in the model appear in this region where 
the chemical potential is sufficiently low and frustrated 
loops are not completely occupied (d < 1). The appear- 
ance of this region is an effect introduced by the disorder 
since in the no frustrated case, it reduces to the point 
fi = 0. 

The entropy per site is 



below which the replica symmetric solution is unstable. 
It can be verified that the above equation is satisfied 
nowhere below Tc and, as in the SK model, here too 
the replica symmetric solution is unstable, although the 
"degree of stability" may vary with fi. 

The susceptibility, x — Pid— q), presents a cusp at 
as can be seen in fig.|. The zero temperature value of y, 
Xo = C /J, depends on /i as shown in the inset of figjj. 
Above ^* it has a constant value, while below it presents 
a maximum. As expected above Tc, x = d/T, depending 
just on the density. We might have chosen to apply h to 
all spins in eq. (|l|) whether their sites were occupied or not 
{hJ2i Si). In this case in the region where d{T = 0) < 1, 
the low temperature susceptibility (now x — /?(! ~~ q)) 
would diverge as due to the free spins which have a 
strong response even to a weak field. 

We report here the values of other quantities in order 
to characterize the system and for comparison with the 
known SG. The internal energy per spin is 



u=\pj'iq^-d') 



J 



d^ 



(19) 



and, at T = 0, 



Uo 



-JCd+^d^ 



djjL 



The compressibility k = P~^dd/dfi has a cusp when the 
transition is continuous, although presents a divergence 
when the transition is first order, as can be seen in figj^. 



s^^P^J\q^-d^) + ^pJd^-p^^d-Pfs . (16) 

AsT^ oo, s 2 In 2 since our phase space has 4^ 
possible states, and at T = 0: 



So = -^C^ + (1 - rf) ln2 



(17) 



For > n* , sq = — l/27r, while for d — 0, sq ^ In 2. Here 
we clearly see the signal of instability of the replica sym- 
metric solution. The next section treats the first step of 
replica symmetric breaking in the Parisi scheme, leading 
to corrections to these results. We also obtain the point 
where sq = 0, that is, /io — —0.026. The reason which 
makes positive the replica symmetric entropy below this 
value is the presence of free spins when the density is 
lower than unity. In those sites where Ui = 0, the spins 
are free to assume any orientation, and this increases the 
entropy (second term in eq.([l7|)). The free energy hes- 
sian eigenvalues have to be calculated in order to verify 
the stability of the replica symmetry solution. Doing so 
we obtain from the condition of positive eigenvalues, the 
AT line M 
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[e~" + cosh{pjyfqz)] 



(18) 



B. Replica Symmetry Breaking 

Results in the previous section have shown the instabil- 
ity of the RS solution and we report here the first Parisi 
correction to it. However, it will appear that the phase 
diagram sketched above in RS is not altered. 

Following the Parisi scheme in the first step of 
replica symmetry breaking (IRSB) the n replicas are di- 
vided in n/m blocks containing m replicas. Different 
replicas in the same block have overlap qi while those 
in different blocks have overlap qo- Thus, the IRSB free 
energy is: 

Pfi = - [(1 - m)ql + mql - d^] - l/JJd^ 

-\n2 - — J Vzo \n J Vzi A™(zo,zi) (20) 

where: 



A{zo, zi) = 1 + e"i cosh(/3/i -I- /3 Jzo\/go + /3JziV9i - Io) 
E, = i^{d-q,) + P{fi-Jd) 

As usual in spin-glass theory, we have to maximize the 
free energy as a function of q, d and m, and the saddle 
point equations are: 
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d=l- Jvzo [A-\ 



qi= J Vza [B% 



(21) 
(22) 
(23) 



(d > qi > qo) and m satisfies 

^m^P^J^ql - ql) - m J Vzo [InAU 

Vzo In / Vzi A™ = (24) 



Here we have defined: 



sinhiPh + (3Jzoy/qo + PJzi^qi - go) 
B(zo,^i) = ^ (25) 



and 



"-^A{zo,zi) 



J Vzi A"^X 



J Vzi A" 



(26) 



The numerical solutions, which can be obtained by ei- 
ther maximizing the free energy or by solving the coupled 
system of saddle point equations, are shown in fig. ^. The 
density does not change in relation to the replica sym- 
metric case, while the effect of replica symmetry breaking 
on qi is rather small, like in the SK model. For m, the 
behavior is analogous to the one found by Parisi in the 
SK model, m = at both T — and T = Tc, and have a 
maximum at an intermediate temperature. 

In analogy with the SK model, the difference between 
fs and /i is only perceptible at low temperatures, al- 
though depending on the value of /i, this effect becomes 
more difhcult to notice. It can also be seen that /i > /s, 
as expected, and that the regions with wrong curvature 
(negative entropy) in the RS solution are largely reduced. 

The susceptibility is given by 



I3[d+ (m - l)qi - mqo] 



(27) 



and the results of IRSB arc higher than the ones found in 
RS (see fig.|), as in the SK model. We have also checked 
that, at T = 0, the point where the density becomes lower 
than one coincides with the one obtained with replica 
symmetry, eq.(p^. 

Our calculations in IRSB confirm that the structure 
of the phase diagram found in replica symmetry is not 
altered when breaking it. The first order transition is 
quite similar to the one found in other spin glass mod- 
els P^-ITtII , although the model deserves a more detailed 
study to clarify whether the transition is or not exactly 
described by one step of replica symmetry breaking. In 
this case the glassy behavior of the model could be rather 
peculiar and different from other disordered system with 
discontinuous transition. 



III. SUMMARY AND CONCLUSIONS 

In brief, we have studied the mean field theory of a 
simple spin glass version of the Blume-Emery-Griffiths 
model with standard replica formalism. This model, de- 
scribed by the Hamiltonian (|^), bridges the usual SG, 
which is obtained in the limit (J//i — > 0, /i — > oo), and a 
version of "frustrated percolation" , which is recovered if 
(^/J^ 0, oo). 

We have seen that the introduction of site variables Ui 
enrich the phase diagram of the spin glass. Specifically 
we have shown that the SG phase around the "frustrated 
percolation" limit, i.e. fj,/J = 0, has peculiarities sig- 
naled, for instance, by a site density below one even in 
the T — limit. Moreover, a discontinuous transition ap- 
pears for negative enough values of /x/J, which may be 
relevant to the description of the structural glass tran- 
sition These two effects, namely the broad d < 1 
region around fi = and the discontinuous transition, 
are not present in the unfrustratcd model where a reen- 
trant phase was present [ p^ . 

The transition is characterized by a cusp in the suscep- 
tibility, although for the compressibility there is a cusp 
when the transition is continuous and a divergence when 
it is first order. By comparing the results obtained with 
IRSB with those found assuming symmetry, we see that 
we are already very near the true mean field behavior of 
the system, since the corrections introduced are small. 
On the other hand, a more detailed analysis should be 
done in the region where the transition is discontinuous 
in order to verify whether, as happens in other models 
with discontinuous transition ||l^-[r^, the IRSB solution 
is exact. 

The mean field results obtained here are in qualitative 
accordance with those found through numerical simula- 
tion for = 3 ||], although a simulation for the infinite 
range version is still lacking and it would be welcomed. 

It would be interesting to study the effect of varying 
the couplings parameters in Hamiltonian eq.(|l|) as 
done in the non frustrated case in [|l^,^ where quali- 
tatively different phase diagrams are so found. Also, by 
introducing a non zero mean in the gaussian distribution, 
it might be possible to study the appearance of the reen- 
trant phase and tune the length of the zero temperature 
< d < 1 interval We expect, by the above results, 
that the interplay of connectivity and frustration would 
lead to a richer static (and dynamic) behavior than the 
usual spin glasses. 
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FIG. 1. The order parameters q and d for several values of 
/u = —0.65 (a), (b) and 3 (c). Notice the diverse behaviour 
of d (and (?) as T — > and the discontinuous transition in (a). 
The respective critical temperatures are Tc ~ 0.13 (a), 0.5 (b) 
and 0.94 (c). 

FIG. 2. The phase diagram T versus /x. The dashed line 

stands for the continuous transition while the solid line corre- 
sponds to the first order one (—0.77 < p < —0.56) and both 
meet at a tricritical point. In particular, when J = — 2/i, 
Tc = J(l + e)"^ ~ 0.27J and when /i = 0, Tc = J/2. 



FIG. 3. The density d versus jj, at T = 0. The point for 
which d-+lis^i = l- (2^)"^/^ while for fj, < -0.77 
there is only the d = solution. 



FIG. 4. The magnetic susceptibility x versus T for = 
and 3 in RS. Inset: T = susceptibility as a function of fj,, 
showing a constant value above /U* and a maximum in the 
region where d{T = 0) < 1. 

FIG. 5. The compressibility /t as a function of n/J for 
T — 0.4, 0.3 and 0.2. In the first two cases the transition 
is continuous and k presents a cusp while in the lEist one it 
has a divergence since the transition is discontinuous. 

FIG. 6. The order parameter qi obtained with IRSB and 
q with RS for = and 3. Regarding d, we did not find any 
difference between the RS and IRSB solutions, while for qi, 
the correction is small. 
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